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Abstract 

Higher dimensional, direct analogues of the usual d = 4 Einstein- Yang-Mills (EYM) systems are 
studied. These consist of the gravitational and Yang-Mills hierarchies in d — 4p dimensional spacetimes, 
both consisting of 2p-form curvature terms only. Regular and black hole solutions are constructed in 
2p + 2 < d < 4p, in which dimensions the total mass-energy is finite, generalising the familiar Bartnik- 
McKinnon solutions in EYM theory for p = 1. In d — 4p, this similarity is complete. In the special case 
of d — 2p + 1, just beyond the finite energy range of d, exact solutions in closed form are found. Finally, 
d = 2p + 1 purely gravitational systems, whose solutions generalise the static d — 3 BTZ solutions, are 
discussed. 



1 Introduction 

Regular and black hole solutions of gravitating gauge field systems have been studied for a long time since 
the pioneering work of Bartnik and McKinnon (BK) PP, where regular solutions of the Yang-Mills (YM) and 
Einstein-Hilbert (EH) systems were presented for d — 4 spacetime dimensions. Other models of gravitating 
non-Abelian gauge fields usually share a number of common features with the BK particles, this example 
becoming canonical. After that discovery, there has been a great deal of numerical and analytical work on 
various aspects of Einstein- Yang-Mills (EYM) theory and a variety of self-gravitating structures with non- 
Abelian fields have been found (for a review see [2]). These include hairy black holes solutions, which led 
to certain revisions of some of the basic concepts of black hole physics based on the uniqueness and no-hair 
theorems. 

In the last years, gravitating solutions with nonabelian fields enjoyed renewed interest in the AdS/CFT 
correspondence, and potentially in the study of Dp-branes of superstring theory. It is therefore pertinent to 
extend the study of such solutions to higher dimensions other than that of d — 4 spacetime. 

Recently regular solutions were found in d — 6, 7, 8 in both regular and black hole solutions in five 
dimensions were presented in 0], the d — 6,7,8 black hole solutions being discussed in [S]. These were 
the classical solutions to systems consisting of higher order terms in both the YM and the gravitational 
curvatures, which do in fact appear in the low energy effective action of string theory. Such terms employed 
in were those constructed from the totally antisymmetrised 2p- forms in both the YM and the Riemann 
curvature 2-forms, namely the YM and the gravitational hierarchies labelled by integers p. The p = 1 
members in each case are the usual YM and the EH systems respectively, while the p = 2 gravitational 
member is the familiar Gauss-Bonnet term. 

These asymptotically flat configurations differ in many respects from the d = 4 BK solutions. In partic- 
ular, all d > 4 solutions have only one node in the gauge potential. Also, the black holes exist only up to a 
maximal value of the event horizon radius. 



The salient property of these solutions [3] 0] is that they exist only for a limited range of the coupling 
parameter a 2 (which gives the strength of the gravitational interaction), and exhibit critical behaviours 
in a 2 . This is not surprising since such composite models necessarily feature more than one dimensional 
constant, analogous to gravitating monopoles Jy. (In the latter case dimensionful Higgs vacuum expectation 
value (VEV) plays this role.) 

More recently, a complete analysis of these critical behaviours in a 2 in the higher dimensional gravitating 
YM systems was presented in |SJ. The methods used in jS] were those of fixed point analysis developed 
previously in j7j for proving the existence of the BK solutions pp to the usual EYM system in d = A. 
Concerning the fairly complicated landscape of critical points here jSj, we restrict our comments only to 
pointing out that in addition to the two types of critical points occurring in gravitating monopoles |6], 
namely the ones associated with the end point and the Reissner-Nordstr0m types, there is also a conical 
fixed point, so described since at that scaling point an angular deficit appears in the gravitational metric 
function. We refer the reader to jS] for a full account of these critical behaviours. 

The analysis in JSj was carried out for models whose gravitational part consisted only of the p — 1 member 
of the gravitational hierarchy, namely the Einstein-Hilbert (EH) term. The physical reason for this was that 
the p > 2 Gauss-Bonnet like terms play only a quantitative role, confirmed numerically in [2], and have no 
effect on the existence of finite mass solutions. By contrast in the gauge field sector, various combinations 
of YM terms with p > 1 essential for the existence of such solutions were employed, as required by scaling 
arguments. Hence higher p > 1 YM terms were employed there, once d was greater than 4. Thus the analysis 
in |S] probed the effect of the higher p > 2 YM terms. 

One of the two main aims of the present work is to probe the effect of higher order p > 2 gravitational 
terms in EYM models. But we know from the numerical results of jS] that once the p = 1 (EH) term is 
present, the effects of all higher p > 2 terms become masked. It is therefore the case that if one wishes to 
study the effects of the p-th gravitational term, all other gravitational terms with pi < p must be excluded. 
In practical terms, this means that we will restrict to models featuring only the p-th gravitational term. 

The other one of our two main aims is to exhibit certain qualitative similarities of the solutions supported 
by a family of EYM models, the first member of which consists of the p — 1 gravitational (EH) and the p = 1 
YM terms, supporting the BK solutions. It follows naturally that the YM term we must choose interacting 
with the (unique) p-th gravitational term, is the (unique) p-th YM term, in spacetime dimensions d = Ap. 
We expect to exhibit a recurring 'symmetry' in the properties of the solutions of this family of models modulo 
Ap dimensions. In the models studied here the complicated critical features of the gravitating monopole [H] 
and of the higher dimensional EYM solutions studied in EJ E] will be absent, since these are due to the 
presence of at least one additional (to the gravitational and YM couplings) dimensionful constant, e.g. the 
Higgs VEV, or, the higher curvature YM coupling constants. 

In Section 2 we introduce the relevant gravitating YM models, and subject them to spherical symmetry. 
These models will be characterised by two equal integers p\ = pi — p specifying the model and the gauge 
group, and the dimension of the spacetime d. In Section 3 we specialise the dimension d to the values 
d = Ap, d = Ap — 1, ...,d = 2p + 2, which are the only dimensions in which asymptotically flat finite energy 
solutions exist. These solutions can be constructed only numerically so the results presented in Section 3 are 
mainly numerical. In Section 4 we specialise to spacetime dimensions d — 2p + 1 in which no asymptotically 
flat, finite energy solutions exist for p\ = p2 = p. The interest in the latter, inspite of their pathological 
properties, is that they can be given in closed form, which is a novel feature in gravitating gauge field theory. 
In addition to these we have supplied two Appendices, A and B, devoted to the extension of some the 
models studied in Sections 3 and 4, to feature a non-vanishing cosmological constant A. Like the solutions 
presented in Section 3, those in Appendices A and B are expressed in closed form. In Appendix A we 
study the geometric properties of the gauge decoupled limit models in d — 2p + 1 appearing in Section 3, 
supplemented with a cosmological term, which are the hierarchy of gravity solutions pertaining to the static 
BTZ solution [£)]. this last being the p = 1 member. The Appendix B presents a generalisation of the exact 
gravity- YM p\ = p2 = p solution in d = 2p + 1 dimensions for a nonzero A. In Section 5, we summarise our 
results. 
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2 The model and imposition of spherical symmetry 



The precise form of the gravitational and non-Abelian matter content of string effective actions beyond 
leading order is still an evolving research subject jl(J|. Our position in studying higher dimensional gravitating 
YM solutions, in 01 Eli has been to choose what we have referred to as the superposition of the p-th 
members of the gravitational and YM hierarchies, consisting of all possible higher order curvature forms 
allowed in any given dimension. 

The said gravitational system is the superposition of all possible (p, g)-Ricci scalars R( P , q ) 

p ^ 

where R( p . q ) are constructed from the 2p-form R(2p) = R A R A ... A R resulting from the totally antisym- 
metrised p-fold products of the Riemann curvature 2-forms R. We express R(p,q) in the notation of as 

j-> _ _/li/i2---Al2p!'ii/2...^_ni n 2 „n q r>m. 1 m2..rn 2p (n\ 

e HP. 1 ?) — vi V2 " Vq t mim 2 ...m2 I ,nin2-.n, J VlP2--W P ' \ ) 

where e" are the Vielbein fields, e = det(e") in 0), and R™^ 2 ''™^ — R{2p) is the p-fold totally antisym- 
metrised product of the Riemann curvature, in component notation. This leads to the definition of the p-th 
Einstein tensor 

G{ P )l = R( P )l ~ ^ e° R {p) . (3) 

One reads from that 

d = 2p + q . 

Now the minimum nontrivial value of q is q — 1, since when q = J5J is manifestly a total divergence, 
namely the Euler-Hirzcbruch density. Thus the highest nontrivial value of P in the superposition Q ls 

Pmax < ~{d- 1). 

The corresponding superposition of the members of the YM hierarchy is 

^-E^^r^p) 2 , (4) 

in which the 2p-form F(2p) is the p-fold antisymmetrised product F(2p) = F /\F A ... A F of the YM curvature 
2-form F. Here the maximum value of P in the superposition Q is simply P max < \d. We define the p-stress 
tensor pertaining to each term in as 

TW = Tr F(2p) AlAlA2 ... A2p _ 1 F(2p)/ lA - A >- 1 - ±- 9fw Tr F(2p) Xl x 2 ...x 2p F(2p) x > . (5) 

The P = 1 systems Q and Q are the usual EH gravity and YM theories, respectively. 

In HUE El > some convenient superpositions (JIJ and Q were selected to be studied, taking into account 
the particular properties of the solutions that were being sought . 

It is our aim here to truncate both |QJ and JIJ such that only one term is present in each. This is so that, 
analogously to the BK case, the solutions not be parameterized by one (e.g. a 2 ) or more parameters. The 
other criterion stated in Section 1 is that in the appropriate dimensions higher members of the gravitational 
hierarchy be employed. So far, subject to respecting the Derrick scaling requirements, one can choose any 
P = pi in and any P = P2 in The final criterion is that of symmetry and analogy with the BK case 
in d = 4, namely that p± = P2 = p in d = Ap dimensions . 

1 There is yet another apparently unrelated coincidence here. In Euclidean signature, the double-self-duality of the Riemann 
2p form curvature leads to the vacuum p-Einstein equations being satisfied, and yields a self-dual SO± (4p) YM 2p-form field 
strength. 
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Thus we define the gravitating YM models in 2p + 2 < d < Ap spacetime dimensions, whose static finite 
energy solutions will be constructed numerically in the next section, by the Lagrangians 

*M = e (§J R (p,o) + 2§fi\ Tr , 2 <1< 2 P> (6) 

k p and t p being constants giving the strength of the gravitational and YM interactions, respectively. Note 
that the gravitational part in © is described by the (p, g)-Ricci scalars R( p , q ) in which q = d— 2p, as defined 
in J2|. For this system, the variational equations for the YM and gravitational fields are 

r p D,(eF^) = 0, « P G W ; = ^T^;. (7) 

The model (01 directly generalises the usual EYM model in d — A spacetime for d = Ap, the latter being 
the p = 1 case. In that case the dimensions of the gravitational part of © are L~ 2p , versus the dimensions 
of the YM part L~ Ap . This choice has been made on grounds of symmetry rather than physics, since the 
leading terms in the effective action of string theory are the p = 1 EH and YM terms, which are both 
excluded. Since the p = 1 EYM system is best known for its BK solutions, the solutions to the hierarchy 
defined by JBJ for d — Ap might be described as the BK hierarchy. 

The definition of the model becomes complete on specifying the gauge groups and their representations. 
Adopting the criterion of employing chiral representations, both for even and for odd d in © , it is convenient 
to choose the gauge group to be SO(d). We shall therefore denote our representation matrices by SO±(d), 
where d = d and d = d — 1 for even and odd d respectively. 

In this unified notation (for odd and even d), the spherically symmetric Ansatz for the SO± (devalued 
gauge fields then reads E] 

.40 = 0, A = (I^>) E<« tl , S <« = 4(I±Ii±l)|I,, rj ]. (8, 

The T's denote the e?-dimensional gamma matrices and 1, j = 1, 2, d — 1 for both cases. 

The spherically symmetric metric Ansatz we use is parameterized by two functions N(r) and er(r) 

ds 2 = -N(r)a 2 (r)dt 2 + N^^dr 2 + r 2 dn'f d _ 2) . (9) 

Inserting (jHJ in ©, the resulting reduced one dimensional Lagrangian is 

f 1 dm \ I 

(10) 



r d - 2 'JP W p 



2-(2p)! 



(2p)N[ -— +(d-2p-l)W 
r dr 



where we note 



(1 - w 2 ) 2 

w= ( - P-. (11) 



3 A Bartnik-McKinnon hierarchy 

It follows from the Derrick-type scaling arguments that static finite energy solutions to the field equations 
of the model © exist only in spacetime dimensions 

2p + 2<d<4p, (12) 

and can only be constructed numerically, which we present here. The d = Ap family is referred to as the BK 
hierarchy, but we study all possible cases allowed by (|12|l . Note that for p = 1, d = A is the only possibility 
allowed by (|P2l . 
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When looking for numerical solutions, it is convenient to redefine 

22 ( P -l) T ^V(4p-2) 



{2 P \)K P 

such that for any (d, p) no free parameter appears in the field equations. 
The field equations implies the relations 



°r r d-2 P -l, 1 _ N yl = r d-2 W p-l 

dr L 



da 

dr 



= 2r 2p - 3 a(l - N) 1 ~ p W p ~ 1 w n 



(13) 
(14) 
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Figure 2. One, two and three nodes solutions of the p = 2 gravity- Yang-Mills model in d — 6 dimensions, 
for the metric function, and 

^-(r^aW^Nw') = ar d - 6 W p - 2 w(w 2 - 1) (^2(p - 1)N(^) 2 + (d - 2p - l)W^j . (15) 
for the gauge potential. 

For globally regular solutions, finite energy requirements and regularity of the metric at r — give 
w = l-br 2 + 0{r 4 ), N= l-4&V+0(r 3 ), a = a (l + 4b 2 r 2 ) + 0(r 4 ), 
where ctq, & are two positive constant. The analysis of the field equations as r — > oo gives 



AT = 1 - 



M, 



i/p 
o 



(7=1 + 



r (d-2p-l)/p 

,2\ PM (1-P)/P 



. .., W = ±l + 



2p{Ac 2 YM y 



r (2p-d+\)/p 
1 



(d-2 P - l)(3p - 1) + 2p 2 r 2p+(d-2p-l)(3p-l)/p 



+ ..., 
+ ..., 



(16) 
(17) 
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Figure 3. The profiles for the metric functions N(r), a(r) and gauge function w(r) are presented for k— node 
globally regular solutions for the p — 3 gravity- Yang-Mills system in d — 12 dimensions. 



with Mo, c arbitrary constants. In analogy to the Einstein gravity, it is natural to identify the constant Mq 
that enters the asymptotic of N with the mass M of the solutions, up to a d— dependent factor. To see this 
let us notice that (fT3| can be written as 

d 



r d-2p-l 



dr 



(1 - NY] 



_ r d-2 T (p)t 



(18) 



i(p)< 



which in the (p — 1, d = 4) case of standard Einstein gravity leads to the usual definition of mass. Here T t 
is the 'energy density' component of the stress-energy tensor defined in (J3J). By keeping with this analogy, 
'the total energy' inside a (d — 2)-sphere Sd-2 with radius r will be given by 

2 



/ dr' / r' d - 2 ^-V) 

'0 JS d -2 



Area(S d -2)r 



d-2p-l 



(1 - N) p (19) 
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Figure 4. One, two and three nodes black holes solutions of the p — 2 gravity- Yang-Mills model in d = 8 dimensions. 



where we made use of the regularity of the function N(r) at origin r = (see (|16H 'I. Taking now the 
limit r — > oo we find that the total energy is finite and given precisely by ~bA=Area{Sd-2)MQ as claimed. 

The equations of motion l|T3 |l -l(TK |l have been solved for p — 2, 3 and a range of d. Considering first 
globally regular configurations, we follow the usual approach and, by using a standard ordinary differential 
equation solver, we evaluate the initial conditions Hlfijl at r = I0~ 4 for global tolerance I0 -12 , adjusting for 
fixed shooting parameter b and integrating towards r — > oo. 

As expected, these solutions have many features in common with the well-known d = 4, p = 1 BK 
solutions. By adjusting the free parameter b appearing in the expansion (|f 6|l . we "shoot" for global solutions 
with the right asymptotics. The solutions are indexed by fc-the number of nodes of the gauge potential w, 
which is always bounded within the strip |u>(r)| < 1. 

For the k-th solution, the function w(r) has k nodes in the interval < r < oo, such that w(oo) = (— f ) . 
Typical profiles for p — 2 solutions are plotted in Figures f , 2 (for d = 8, 6) and p = 3, d = 1 2 configurations 
(Figure 3). It is worth noting here that for a p-gravity and p-YM theory, when the spacetime dimension is 
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Ap, then the profiles of the functions N(r), a(r) and w(r) are perfectly smooth as seen from Figures 1 (for 
d = 8) and Figures 3 (for d = 12), as is the case for the BK solutions in d — A, while for dimensions different 
from d = Ap, e.g. for d = 6 depicted in Figure 2, these profiles are less smooth and the numerical analysis is 
correspondingly more delicate. This is one the features highlighting the similarity of the solutions to these 
theories modulo 4p. 

The behaviour of the metric functions N and a is similar for all fc's. The metric functions a increases 
with growing r from tr(0) = do > at the origin to <t(oo) = 1. As k increases, N develops a more and more 
deep minimum at some r m , closely approaching the zero value for large enough values of k, as indicated in 
Figures 1-3. At the same time, the value at the origin of the metric function a strongly decreases with k. 

The limiting solution with k — oo can also be investigated |12| . Since the Schwarzschild coordinate system 
breaks down in this limit, we should use a different parameterization and the equations are formulated as a 
dynamical system. Similar to the p = 1 case, this turns out to be non-asymptotically flat. 

The system <|13[) - (I15[) presents also black hole solutions. The boundary conditions at infinity are still 
given by l|17[) , and are now supplemented with the requirement that there is a regular event horizon at some 
r = rh > 0, with N(r) > for r > r^. The local power series in the vicinity of the horizon reads 

w(r) = w h + w' h {r - r h ) + 0(r - r h f , N(r) = N' h {r - r h ) + 0{r - r h ) 2 , a{r) = a h + o' h {r - r h ) + 0{r - r h ) 2 , (20) 

with 

N' h = {d - 2p - l \ l-r^Wl), w > h = Wh (wi-l)(d-2p-l)^-, a h = 2a h rf-*W»- 1 w'l (21) 

/-^ ^2 \ 2 

and Wh — ~4 ■ Similar to the d = 4, p = 1 case, one find a sequence of global solutions in the interval 
r^ < r < oo, for any value of the event horizon r/> > 0. These solutions are parameterized again by the node 
number k of the gauge function w. For any (r^, k) the behaviour of the functions w, a is qualitatively similar 
to that for regular solutions. The gauge function w starts from some value < Wh < 1 at the horizon and 
after k oscillations around zero tends asymptotically to (— l) fe . In this case again one can show that \w\ < 1 
everywhere outside the horizon. In the limit — > the event horizon shrinks to zero and the black hole 
solutions converges pointwise to the corresponding regular configuration. In Figure 4 we exhibit the one, 
two and three nodes solutions of the p = 2 model in eight dimensions; similar solutions have been found for 
p = 2, d = 6 and p = 3, d = 9, 12. 

The stability of solutions in the usual p — 1, d = 4 case has been studied extensively, both perturbatively 
and at the non- linear level (see the discussion in and the references therein). It turns out that all 
known regular and black hole solutions in that case are unstable with respect to small spherically symmetric 
perturbations. The most obvious indication of the instability comes from the absence of a topological charge 
in the YM sector. This is obvious in the d = 4 case, but when several members of the YM hierarchy are 
present like in H3EJIH1, in some of these theories a Pontryagin charge is defined, leading to stable solutions. 
The stability question in such models has been studied in J3| • I n t ne present work however we have restricted 
to a single (the p-th) member of the YM hierarchy, so that no topological charge can be accommodated in 
the case of static solutions. Indeed, the stability analysis in J2j leading to the conclusion that the solutions 
in dimensions in which a Pontryagin charge is not defined are unstable (sphalerons) applies to the models 
studied in the present Section, and is not repeated here. 

Another aspect of the stability of gravitating YM fields is the interesting case when a negative cosmological 
constant is introduced, in the usual EYM model with p = 1, d = 4 |141 115j . In that case, one can see from 
the asymptotic analysis that the value of the gauge field function w(r) at infinity is not restricted to ±1, 
and in particular linv^oo w(r) = is allowed. This is the asymptotic value for a static monopole, i.e. 
the particular solution in question is a finite energy lump with (topological) monopole charge, rendering 
it stable. Unfortunately this property does not persist 2 in the p > 2 models with negative cosmological 
constant, implying that these would be exclusively sphalerons. As a result we have eschewed an analysis of 
these. 

2 This can be seen from the large r asymptotic analysis in the case of p > 2 models with negative cosmological constant, 
which leads exclusively to lim 7 ._ (00 w(r) = ±1. 
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We close this Section with several comments on the thermodynamic properties of the black hole solutions. 
The Hawking temperature of these configurations can easily be found by using the standard Euclidean 
method. For the line element ©, if we treat t as complex, then its imaginary part is a coordinate for a 
non-singular Euclidean submanifold iff it is periodic with period 

47T 

= N'(r h )<r(r h )' (22) 

Then continuous Euclidean Green functions must have this period, so by standard arguments the Hawking 
temperature is T = 1//3. For any (p, d) this is found to decrease with the node number k. To compute the 
entropy of these solutions, one may use the standard relationship between the temperature and the entropy, 
S = J dM/T + Sq, with So a mass- independent constant. Although further study is necessary, we expect 
this expression to differ from the standard one quarter event horizon area value, which holds for p = 1 EH 
gravity only ^| E| • 



4 Exact gravity- Yang-Mills solutions in d = 2p + 1 

In the particular dimensions d = 2p + 1 , the analogues of © , 



C 



(2p+l) 



R 



2p " (p ' 1} ' 2(2p)! 

do not support static finite energy solutions, however, their solutions can be constructed in closed form. 
Since to the best of our knowledge no exact (nontrivial) solution is known in the literature for the coupled 
gravity- Yang-Mills equations 3 , we discuss here the basic properties of these d = 2p+ 1 configurations. The 
straightforward generalisation in the presence of a cosmological constant is presented in Appendix B. 
After several redefinitions of the theory's constants, the effective Lagrangean of this theory reads 

2 - 



TrF(2p) : 



(23) 



L = a 



d 
dr 



(N - l) p - cr^Niw 2 



1 



dw 
dr 



(24) 



where c ~ t/k is a free parameter (in this general case, we do not fix the sign of c). It is also convenient to 
redefine the gauge potential according to 



a(w) = (w 



i)p-i dw = (-ir +i 2 F 1 (±i-p,^w 2 ), 

(2-Fi(a, 6, c, z) being the hypergeometric function), such that the system 124fl admits the first integral 



{w 2 



1) P ~V 



,l-2 PjVcr ' 



where a is an arbitrary real constant. 

The metric variables a and N satisfy the equations 

dXP 



dr 

dY 

dr 



ca 



(X + 1)Y 



2ca 2 



r 2p-l 



p {x + ifxp- 11 



where X = N — 1, Y = a 2 . This implies the relation: 



dY 
dX 



2Y 



(X+1)' 



which gives 



Y(X + l) 2 = C, 



(25) 
(26) 

(27) 
(28) 

(29) 
(30) 



3 See, however, the exact solution with planar symmetry of the d = 4 EYM equations with a negative cosmological constant 
A = —3 presented in 1181 . 
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Figure 5. Penrose diagram of the metric HM6^ for c > 0. Wavy lines correspond to curvature singularities. 

(i.e. No = C 1 / 2 ) and we can set C = 1, without any loss of generality. 4 Replacing this expression 
into the X-equation one finds the relation: 



X P 2 F 1 (p,l 1 l+p,~X) = — 

2p 



,2p 



■P, 



(31) 



where (3 is an integration constant. 

Expressed in a form which employs X as coordinate, the general solution reads 



where 



ds 2 = h(X)dX 2 



X + l 



pXP~ 



ca 2 (X + 1) 



r 2-4p 



r(X) 



^( X P 2 F 1 (p,lA+P,-X)-P) 



It follows straightforwardly that the expression for the transformed gauge potential a is 

a(X) = a a + -r 2 P(X), 
2p 



(32) 

(33) 
(34) 

(35) 



with do a constant of integration, which from (|25|l gives the expression of the gauge potential. One can see 
that this general solution is not asymptotically flat (as the matter fields do not vanish at spacelike infinity). 
The solution expressed in the r-coordinate takes a simple expression for p = 1 only 



ds 2 



foe' 



\dr 2 -dt 2 ) + r 2 d(p 2 



w(r) — w + 



(36) 



where a, /o and wq are integration constants. In order to eliminate the conical singularities in the (r, ip) 
sector we require that fa = 1. It is straightforward to compute the Kretschmann scalar for this metric: 



A' 



1 ^^iva p 1 L 



3c 2 a 4 e c 



/o 2 



(37) 



'We assume that the integration constant C is positive, i.e. a 2 > 0. 
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while the Ricci scalar is proportional to e ca r / 2 . Notice that if c > these curvature scalars blow up in the 
limit r — > oo, while apparently the geometry is well behaved and becomes regular when c < 0. The Penrose 
diagram of the metric H'UiH in the c > case is presented in Figure 5. Each point in the diagram represents 
a circle of radius r, as described by the ip coordinate. 

While it might be surprising to find null curvature singularities at infinity, to understand their presence 
notice that defining the following null coordinates: 



we can write the metric as: 



while K becomes: 



u = t — r, v = t + r, (38) 



ds 2 = -f e- ca2{v ~ u)2/8 dudv + ^(v~u) 2 d<p 2 (39) 



Q r 2 4 ca 2 0-u) 2 /4 

K = (40) 

Jo 

It is now clear that in the (r, t)-plane null rays are represented as straight lines at 45° with respect to the 
coordinate lines and furthermore K diverges as u — > oo with v = const, as well as for v — > oo with u = const.. 

To get a better appreciation of the properties of the c > 3-dimensional geometry it is of interest to 
examine in more detail the behaviour of the timelike and null geodesies. These will illuminate the nature 
and the effects of the naked curvature singularities at infinity. Since t and tp are cyclical coordinates we can 
write down directly the following constants of motion: 



2 2 , 

g — ca r ' 



/2 i = E, r 2 if = L (41) 
while the radial coordinate r satisfies the equation: 



,2 



-r 



V(r) = (42) 



which is the equation of a material point with unit mass and zero energy, moving in an effective potential 
given by: 

£V qV LV qV / 2 ee caV / 2 , , 

V ^ = 2— + ^ 2— (43) 

where e = — 1 for timelike geodesies, e = for null geodesies and e = 1 for spacelike geodesies. Consider 
first the timelike geodesies, i.e. e = —1. For non-radial geodesies L ^ and for E ^ we have V(r) — > oo 
and the region near origin is classically forbidden, while in the limit r — > oo we have V(r) — > — oo, which 
means that the particles accelerate so that \r | — > oo as r — * oo. For radial geodesies L = and the potential 
reaches a finite value as r — > 0, while again V(r) — > — oo as r — > oo. These properties continue to hold for 
null geodesies. Let us notice that null and even radial timelike geodesies can reach spatial infinity in finite 
intervals of the affine parameter. Thence the c > geometry is pathological. 

By contrast, if c < 0, generically V(r) — > when r — > oo and therefore the region r —> oo is indeed at 
infinity in terms of proper distance and null right rays cannot reach spatial infinity in finite intervals of the 
affine parameter. Then for c < the geometry is indeed free of any curvature singularities. That this is 
indeed the case can also be confirmed by examining the components of the Ricmann tensor components in 
an orthonormal frame in the large r limit. 

For either sign of c, near origin the spacetime geometry is flat and regular if /o = 1. It is interesting 
to note that one can also have the case in which /o < 0, in which case the geometry is time-dependent. 
However, one can easily check that the geometry is pathological if c > and regular for c < 0. 
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One should also notice that the 3-dimensional metric is directly written in a Weyl-Papapetrou form. 
When c < this will allow us to make an interesting connection with a specific time-dependent axially 
symmetric geometry of the form: 



ds\ 



e 



\/-2ea 2 i j 2 



dz + e 



V-2ca 2 t 



foe 



/2 (dr 2 -dt 2 )+r 2 d(p 2 , 



(44) 



which is a solution of the vacuum EH field equations in 4 dimensions, = 0. Upon dimensional reduction 
along the z coordinate we obtain precisely the 3 dimensional geometry (|36|l that is now supported by a 
time-dependent scalar field <fr = v 7 — 2ca 2 t, which is a solution of the equations of motion derived from the 
Lagrangian: 



One can now perform the dualisation of the scalar field to obtain a magnetic 2-form field strength and it is 
easy to check that the final solution corresponds to the magnetic version of the Reissner-Nordstr0m solution 
in 3 dimensions |19 | -|21 | . 

The 4 dimensional geometry is also free of curvature singularities. If we perform the analytical continu- 
ations t — > iz, z — ► it and also c — * — c we obtain the Euclidean version of the 3 geometry (|36|l where now 
c > 0. Its Lorentzian version is obtained if we further analytically continue z — > it with the final result that 
JHEJ with c> is now supported by a time dependent scalar field whose kinetic term has the wrong sign. 
This can be regarded as an indirect confirmation of our previous result that the geometry il36|) for c > has 
pathological properties. 

5 Conclusions 

The general aim of this work was to study the properties of gravitating gauge field systems whose gravitational 
part consists of higher order gravitational curvature terms, e.g. Gauss-Bonnet terms, i.e. necessarily in higher 
dimensions. In particular, our scope is limited to static finite energy solutions, which are also spherically 
symmetric and asymptotically flat. As we know from previous work that the presence of the usual Einstein- 
Hilbert (EH, p = 1 gravity) term masks the effects of all higher order (p > 2) terms, we have chosen to 
study models featuring only one gravitational term with p > 2, whose specific value is chosen according to 
the dimensionality d of the spacetime. This still leaves open the choice of the YM terms. 

A dual aim in this work was to bring out some general features of gravitating YM solutions, seen in the 
original BK pQ case with d = 4, EH p = 1 gravity, and usual p = 1 YM term. This has led us to study models 
with p-th order gravity and 2p-th order YM in d — Ap dimensions. This is the content of results in Section 3, 
where we have verified that all qualitative properties of the BK solutions are repeated in dimensions modulo 
Ap (although nontrivial finite mass solutions exists also for 2p + 2 < d < Ap). Like the BK solution, these 
are constructed numerically and are likewise all (unstable) sphalerons. 

As a result of our general considerations we realised that a somewhat different hierarchy of models, namely 
those again with p-th order gravity and 2p-th order YM but now in d — 2p+ 1 dimensions, actually supported 
exact solutions in closed form. These configurations were discussed in Section 4. Their properties depend 
essentially on the sign of the coupling constant c ~ r/k. While we expect that the physically meaningful 
configurations be described by negative values for c, in Section 4 we studied in more details the p — 1 case. 
We found that indeed c > leads to a pathological geometry, while for c < the solution is perfectly regular 
and we anticipate these properties to hold for higher values of p. Led on by the work in Section 4, featuring 
solutions in closed form, we branched out to introduce a cosmological constant to the models studied in 
Section 4. This also led to solutions in closed form, presented in the Appendices A and B. In Appendix A we 
considered the gauge decoupled versions of these models since this turns out to yield a hierarchy of solutions 
the first (p = 1) member of which is the static BTZ solution Unfortunately the remarkable geometric 
features of the p = 1 solution [Hj are not repeated in the p > 2 cases, and a fairly extensive analysis of this 
is given in Appendix A. We found that for generic values of the parameters these solutions are pathological 



£ 3 = ei?, 



'■(1,1) " \ e i. d 4>) 2 



(45) 
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in that they exhibit naked curvature singularities (that can be hidden inside cosmological type horizons), 
however, for special values of the parameters regular (a)dS '/ flat backgrounds are obtained. In Appendix 
B, the closed form solutions of the models in Section 4 with a non- vanishing A are presented. The general 
solution has two branches and in the A — > limit only one of them will survive to give the solution discussed 
in Section 4. The general form of these solutions is very complicated and this impedes a general analysis of 
their properties. We do expect however that since they have no horizons, they will generically exhibit naked 
curvature singularities. 

While one might question the physical relevance of the new exact solutions found in this paper since the 
form of our gravitational and matter Lagrangians is non-standard, we take the point of view that given the 
scarcity of known non-trivial exact solutions of the YM system coupled to gravity, any new exact solutions 
that can be found in closed analytical form might shed some light on the properties of such complicated sys- 
tems. Moreover, in literature there have been studied modifications of the standard Einstein-Hilbert gravity 
by considering invariant quantities constructed from the curvature scalar and/or the Riemann tensor. In 
general it is known that for a gravitational Lagrangian constructed out of the metric and the Ricci tensor it 
is possible to perform a conformal transformation (or more generally a Legendre transformation) to a metric 
expressed in the Einstein frame, solution of the standard Einstein-Hilbert gravity coupled with (exotic) mat- 
ter fields (see for instance 22 23] and references therein). The question if a similar reasoning can be applied 
to the gravitational lagrangians considered in this paper remains an interesting topic for further research. 
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A A BTZ hierarchy in (2p+ l)-dimensions 

It is well known that while gravity in d = 3 spacetime dimensions in the absence of matter is dynamically 
trivial, in the presence of a cosmological constant it supports highly nontrivial asymptotically anti-deSitter 
(AdS) solutions (BTZ) [§]. 

Our purpose here is to present a hierarchy of gravitational models with non-vanishing cosmological 
constant, which support exact solutions in closed form, proposing a generalisation of the static d = 3 BTZ 
solution. It turns out that these models are defined in spacetime dimensions d = 2p + 1 and are described 
by © with the YM terms suppressed. This can be written by augmenting ©, with YM fields suppressed, 
by a cosmological constant, 



(2p+l)!eA, 



(A.1) 



-(2p+l)| F=0 

whose p = 1 member supports the familiar BTZ solution [[]]. 

The exact solutions in question are the static spherically symmetric field configurations resulting from 
the imposition of spherical symmetry, given by the metric Ansatz in Schwarzschild coordinates, where 
here r is the 2p dimensional radial coordinate. 

Here, we introduce the cosmological constant, i.e. we adopt the system (|A.1|) . Subjecting the latter to 
the Ansatz ©, after a suitable rescaling we have the reduced Lagrangian 5 

(d-2)\ 



1 



2 2 p~! (d-2p 
We immediately find the following solutions 

1 



1)! 



dr 



(N - If - r 2p - x A 



{N -If 



2jg_ 



r 2p A + const. 



p{N- \y 



_i da 
dr 







=> a = const. 



(A.2) 



(A.3) 



5 Note that A has here the opposite sign as compared to the standard choice in literature 
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Notice that for even values of p we have two branches in our solutions. We will discuss the global structure 
of these metrics separately according to even or odd values of p. Also, we will denote the constant that 
appears in (|A.3(1 by M, as this integration constant will be proportional in some cases to the mass of a black 
hole for black hole type spacetimes. We will also define A = ±^ to simplify notations. For A > we obtain 
then in general: 

N(r) = 1±I^— + M) , (A.4) 

where the minus sign defines a second branch of solutions for even values of p only. To analyse the case of 
a negative cosmological constant A < one simply analytically continues I — > i£, while the case A = is 
obtained in the limit £ — ► oo. 

Several particular cases of interest are p — 1 and p — 2. The former solution corresponds to the celebrated 
BTZ solution for A = Jr > 0: 

2 



ds 2 = -N(r)dt 2 + N-'{r)dr z + r 2 d8 



r 2 



N(r) = — — (A.5) 

while A < corresponds to 3-dimensional dS space. If A = the spacetime is flat and, unless M = 0, it 
contains a conical singularity at origin. For p = 2 the metric is 5-dimensional and solves the pure Gauss- 
Bonnet equations with a cosmological constant: 



ds 2 = -N{r)dt 2 + N~ 1 (r)dr 2 + r 2 dQ 



3 



f 4 



N(r) = l±\±j i +M) (A.6) 

In general, we have two branches of our solutions, which correspond to the choice of the sign in front of 
the radical. Inside the radical the choice of sign is dictated by the sign of the cosmological constant. For 
a negative cosmological constant we have to consider strictly positive values for M and, in this case, the 
radial coordinate r will take values only in a finite interval. We find, however, that there are curvature 
singularities at both end points of this interval, separated by a black hole type horizon in between. For a 
positive cosmological constant, both positive and negative values of the parameter M are allowed as long 
as the expression under the radical is positive. We find that there is a naked curvature singularity located 
at r = 0, which is hidden inside a cosmological type horizon. When M — we obtain AdS spacetime as 
the positive branch solution, respectively dS as the negative branch solution. If A = we obtain a solution 
containing a naked curvature singularity at origin and having a deficit of solid angle. As we shall see in 
the followings, the general solutions in higher dimensions exhibit similar properties as the lower-dimensional 
cases. 

The detailed analysis of the global structure of these spacetimes involves a discussion of the singularities, 
horizons and the asymptotic structure. We will be mainly interested in curvature singularities and these 
will be identified using the Kretschmann invariant. For the spherical symmetric ansatz (j5J) the Kretschman 
scalar can be written as: 

K — R rwp 

= (Ny + 2 -^(N>) 2 + 2id -^ d - 2 \l-N) 2 (A.7) 

From the general form of N(r) we can see that there is a curvature singularity located at r = and also 
at points where A + M — 0. Around these curvature singularities the Kretschmann invariant behaves as 

K ~ O ( ) . Notice that these metrics are regular only if M = 0. 
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The asymptotic behaviour is controlled by the dominant term in N(r) as r — > oo. In general, for odd 
values of p one can have asymptotically (a)dS spaces, however, for even values of p we find that we have to 

2p 

restrict the range of the radial coordinate such that ±?3p + M > 0. For A = we find that the asymptotic 
structure is controlled by the values of the parameter M. If M = we obtain the flat spacetime. If M ^ 
then in general we obtain spaces with deficits or surfeits of solid angle. 

To characterise the horizons we use the following definitions |24| : a horizon located at r = r% is a null 
hypersurface with finite curvature, such that N(rh) = 0. A black hole horizon is defined by the condition 
N'{rh) > 0; a horizon for which N'(rh) < and r/, is the largest root of N(r) will be called a cosmological 
horizon. If N'(rh) < and r/i is not the largest root then r = 17, will define an inner horizon. If N'(rh) = 
then r = would correspond to an extreme horizon. Notice however that in our solutions this can happen 
only if rh = and since r = is the location of a curvature singularity the spacetime will be singular. 

Consider first spacetime geometries corresponding to odd values of p. We will examine the global structure 
for each value of A and M. 

• A > For any value of M the spacetimes are asymptotically adS. If M > 0, there is a curvature 
singularity at r = and since N(r) does not vanish for any value of r there are no horizons, hence 
the spacetime contains a globally naked singularity. If M — then the spacetime is adS. For M < 

r 2p 2 p 

we find that N(r) can have a zero at r — ru such that -fa + M = — 1. However, since + M = 
for a value r s > rn, then r — r s defines a curvature singularity that is not covered by a horizon and 
therefore it corresponds to a globally naked singularity. 

• A = In this case N(r) = 1 + Mp and there is a curvature singularity at r = 0. We exclude from 
our discussion the value M = — 1 for which N(r) = 0. If M > by rescaling the coordinates we can 
bring the metric in the following form: 

ds 2 = -dt 2 + dr 2 + (1 + Mp)r 2 dQ 2 d _ 2 (A.8) 

The spacetime has a surfeit of solid angle as one can see by computing the surface area of a sphere 

with radius r and comparing its value with the one calculated in flat spacetime. If M = we obtain 

i 

the flat spacetime. For negative values of M such that 1 + M 5 > the spacetime has a deficit of solid 
angle. If 1 + Mp < then r is a timelike coordinate and the spacetime corresponds to a Milne-type 
spacetime having a deficit (or surfeit) of solid angle for M < —2 P (respectively M > — 2 P ). 

• A < If M < then N(r) — when ™£ — M = 1 and the curvature singularity at r = is hidden 

behind a horizon located at ru = 1(1 + M)>. Notice that in order to obtain real values for we have 
to restrict the values of M such that M > — 1. Since N'(rh) < this corresponds to a cosmological 
horizon and the spacetime is asymptotically dS. For M = the spacetime becomes dS while for M > 
the spacetime contains a globally naked singularity located at r s — lM*p . 

Let us consider next the geometries corresponding to even values of p. As we have previously mentioned, 
there are two branches of solutions and we shall discuss each branch separately. The negative branch has: 

N(r) = l-i±— +M j (A.9) 

Here the upper sign corresponds to a positive cosmological constant A > 0, the lower sign corresponds to 
A < 0, while the case A = is obtained in the limit I — > oo. Notice that we have to restrict the values of 

2p 

the parameters M, I and of the radial coordinate r such that ±725- + M > 0. 

• A > For M > the spacetime is asymptotically dS and it has a curvature singularity at r = 0, 
hidden inside a cosmological horizon located at = 1(1 — M) *p . Notice that in order to obtain real 
values for r/, we must restrict the values for M such that M < 1. For M = 1 then r/, = and the 
spacetime contains a globally naked singularity. If M > 1 there is a globally naked singularity located 
at r = Q. If M = we obtain dS. For M < we find a singularity at r — inside of a cosmological 
horizon located at — l(\ — M)*p . 
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• A = In order to avoid complex values for N(r) = 1 — M> we must require that M > 0. If M = 
we obtain Minkowski spacetime. If M > the spacetime contains a globally naked singularity at r = 
and it has a deficit of solid angle. 

• A < If M < we find that /(r) takes complex values. If M > we must also restrict the values of 
the radial coordinate to a finite interval such that M — > 0. There exists a black hole type horizon 

at Th — i(M — 1)5F only if M > 1. However, at both endpoints of the radial coordinate there are 
curvature singularities. For M = 1 there is a naked curvature singularity at r — and also another 
one at r = £. 

Consider next the positive branch solutions. In this case we have: 

N(r) = 1+U—+M) (A.10) 
The analysis is very similar with the one performed for odd values of p. There exists however an extra 

2p 

restriction that ±jjf + M > 0. The equality sign corresponds to a curvature singularity location. We find 
then that these spacetimes are generically singular unless M — and in the cases where there are regular 
horizons, these correspond to cosmological horizons such that the spacetimes still contain globally naked 
singularities. 



B D = 2p + 1 gravity- Yang-Mills solutions with A ^ 



For p = 1, d = 4, the EYM equations with a cosmological term present solutions with very different properties 
as compared to the A = case ^1 ■ Their d > 4 generalisations in Einstein gravity with higher terms in 
the Yang-Mills hierarchy have been discussed recently in It is therefore natural to consider solutions of 
the p-th. gravity- Yang-Mills system l|23|l in the presence of a cosmological constant, in d = 2p+ 1 dimensions. 
Using the same notations as in Section 4, the reduced Lagrangean of this system reads 



L = a 



(N - If - cr^Niw 2 - ifP- 2 (j^j ' 



(B.l) 



the resulting field equations being solved by using the same techniques as in the A = case. 

Employing the same variables, X = N — 1 and Y = a 2 , we find the same c-equation, while the function 
N satisfy the equation 



dX p 
dr 



r 2p-l 



ca 



{X + l)Y 



A 



(B.2) 



the gauge field equations still admitting the first integral l|26ll . 
Therefore Q29J1 is generalized to 



dY 
dX 



2Y 



(X + 1)(1 + AF(X + 1))' 
where A = A/ca 2 . If we define a new function Z = (X + 1)Y, we find that 1B.3|) can be written as: 

Z XZ - 1 



and upon integration we obtain: 



dZ 
dX 



(XZ-1) 



X + 1 XZ + 1 



X + l 



c 



(B.3) 



(B.4) 



(B.5) 
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where C is an integration constant. This equation has now two solutions, denoted by Z± and given by: 

\Z± 



\/x + ci =f VaTi 



y/X + Ci ± VX+T 
where c\ = 4AC + 1. In consequence, we obtain two branches for our solution 

• a vx + ci ± VaTI 



ccr 



± 



A(l + X) yx + ci T VXTT' 



(B.6) 



(B.7) 



Notice that when taking the limit A — > only the negative branch will survive and we obtain precisely 
the solution given in (|30|) . The positive branch solution has no asymptotically flat limit. Replacing these 
expressions in Q27[l. we arrive at a relation similar to (|31|l : 



,2p 



2p 



where 



F±{X)= J XP- 1 U± 



X + l 
x + Cl 



dx. 



(B. 



(B.9) 



It appears that it is not possible to find a general expression of this integral. Several particular cases of 
potential interest in which we can integrate HB.9|I are: 



F±(X) = X ± y/(X + TyX + ci) T (ci - 1) log(VxTT + ^A + ci), for p = 
F±(A) = -j-(4X 2 ± 2y>{X + 1)(X + ci)(l - 3ci + 2X) 



1. 



(B.10) 



±(ci - l)(3ci + 1) log(l + ci + 2X + 2y/(X + l)(X + ci))) for p 
The general solution takes a simpler form when expressed using X as coordinate with 



r(X) = \^ K F ± (X)-2pp 
Therefore, the general metric of the A ^ solution is given by 



with 



ds 2 = 

9i(X) = 
92(X) = 



gi {X)dX 2 + r 2 (X)dft 2 l _ 2 - g 2 (X)dt 2 



pr -2 P -i (x) 
2A 



1± 



X + l 
X + Cl 



X 



(B.ll) 



(B.12) 



(B.13) 



za 2 + ci ± y/X + 1 

Ty/x + dT VxTT' 



From the above form of the general solution, notice that we have to restrict the X coordinate such that 
X > —1 and that there are no horizons. However g\{X) will blow up as X — > —1 and therefore we conclude 
that in general such spaces will have pathological features. 
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